Magnetic translation algebra with or without magnetic field 
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The magnetic translation algebra plays an important role in the quantum Hall effect. Murthy 
and Shankar have shown how to realize this algebra using fermionic bilinears defined on a two- 
dimensional square lattice. We show that, in any dimension d, it is always possible to close the 
magnetic translation algebra using fermionic bilinears, be it in the continuum or on the lattice. We 
also show that these generators are complete in even, but not odd, dimensions, in the sense that 
any fermionic Hamiltonian in even dimensions that conserves particle number can be represented 
in terms of the generators of this algebra, whether or not time-reversal symmetry is broken. As 
an example, we reproduce the /-sum rule of interacting electrons at vanishing magnetic field using 
this representation. We also show that interactions can significantly change the bare band width 
of lattice Hamiltonians when represented in terms of the generators of the magnetic translation 
algebra. 



I. INTRODUCTION 

The two coordinates of an electron in the plane orthog- 
onal to an uniform magnetic field B e pointing along the 
direction e do not commute. There follows the U(l) al- 
gebra 



[t(a),t(b)} 



(a A b) 



t(a + b) (1.1) 



obeyed by the triplet of generators t(a), t(b), and t(a+b) 
for any pair a and b of vectors orthogonal to e, which is. 
called the magnetic translation algebra in this context .□ 
The magnetic translation algebra can be used to derive 
the transverse conductivity of the integer quantum Hall 
effect (IQHE). It has also been used by Girvin, Mac- 
Donald, and Platzman in Rcf. |^ to place a variational 
estimate on the excitation gap for the fractional quan- 
tum Hall effect (FQHE), following closely the approach, 
of Feynman and Bijl in their study of excitations in 4 He.u 

Hamiltonians defined on two-dimensional lattices with 
topologically non-trivial bands can also display quan- 
tum Hall physics. The IQHE can occur in band insu- 
lators when the Bloch bands haw a nonvanishing Chern 
number, as shown by Haldane.Q The FQHE effect re- 
quires strong electronic correlations. This is possible, if 
the Chern bands are sufficiently narrow (or even flat).ErEI 
Whether flat Chern bands can sustain or not a FQHE is 
a matter of energetics. Exact diagonalization studies of 
fractionally filled Chern bands with added short-range in- 
teractions are consistent with a correlated liquid ground 
state supporting a FQHE for certain filling fractions. ErE2l 
Such topological correlated states on the lattice are now 
known as fractional Chern insulators (FCI). 

In an effort to draw a bridge between the case when 
the FQHE is realized in the continuum or in a FCI, 
Parameswaran, Roy, and Sondhi in Ref. ^l] have pio- 
neered an algebraic approach to FCIs by deriving the 
algebra obeyed by the-detisity operators projected to the 
partially filled band.U'Eacil They found that the alge- 



bra (Id) emerges to leading order in a gradient expan- 
sion. Remarkably, Murthy and Shankar have (i) con- 
structed in Ref. |25| a coherent superposition of the pro- 
jected density operator that closes the U(l) algebra < \l . l| ) 
on the square lattice and (ii) represented any Hamilto- 
nian that commutes with the number operator and de- 
scribes the competition between the electronic hopping 
and the electronic interaction in terms of these genera- 
tors. 

In this paper, we are going to generalize the results 
by Murthy and Shankar as follows. We shall represent 
the U(l) algebra QL.lD in terms of coherent superposi- 
tions of electron-hole pairs in arbitrary dimensions both 
in the continuum and for Bravais lattices. We shall then 
show that these generators provide a complete basis for 
the linear space of operators spanned by charge neutral 
fcrmion bilinears provided the Bravais lattice, or its em- 
bedding space in the continuum limit, is even dimen- 
sional. For odd dimensions, the generators of the U(l) 
algebra (1.1) form an incomplete basis of the space of 



operators spanned by charge neutral fermion bilinears. 

We treat the case of Hamiltonians acting on wavefunc- 
tions supported in the continuum in Sec. O. We then turn 
our attention to Hamiltonians actin g on wavefunctions 
supported on Bravais lattices in Sec. III. We verify that 
the /-sum rule is obeyed if one represents the electronic 
density operator in terms the particle-hole generators of 
the algebra (1.1) in Sec. IV A. In Sec. IV B, we explain 
how, when represented in terms of these generators of the 
algebra (1.1), interactions induce one-body terms that 
can significantly change the bare band width of lattice 
Hamiltonians. The same effect in the FQHE requires the 
addition of a strong one-body perturbation to a Landau 
band, one that is of the order of the FQHE gap. Thus, 
whereas the FQH is a strong coupling problem, the FCI 
in a flat band is more like a problem at intermediate cou- 
pling. This result explains why in Ref. [t8] an FCI with 
Chern number two was more stable if the bare dispersion 
was not flat rather than flat, for the bare and induced 
one-body terms can conspire to neutralize each other. 
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II. 



THE CASE OF THE CONTINUUM 



that obeys the algebra 



We define the fermionic Fock space $ with the help of 
the algebra 



{c(k),c^(k')} = S(k-k'), 
{c(k),c(k')} = {c\k),c*(k')}=0, 



(2.1) 



for any pair of momenta k, k' <G M. d . Without loss of 
generality, we ignore any internal degrees of freedom such 
as the spin quantum numbers since we are after the U(l) 
algebra (1.1). 

The linear space of fermionic bilinears that we study 
is spanned by the basis 



T( qi ,q 2 ) :=c\ qi )c(q 2 ) 



(2.2a) 



f(q 1 ,q 2 ),f(q' 1 ,q' 2 )\=5(q 2 -q' 1 )f(q 1 ,q 2 ) 

-S{q 1 -q! l )f(g' 1 ,q 2 ) 



(2.2b) 



for any quadruple q 1 , q 2 , q[ , and q 2 from R d . 

For any momentum q £ K d and for any function / : 
R d x R d — > C, define the coherent superposition 

Q f (q)~ J f(q,P)c t (q + P )c(p). (2.3a) 
p 

There follows the algebra 



'Q f (q),Q f '(q') 



f(q, q' + p) /'(«', P) ~(q^q' and / O /') c\q + q'+p)c (p) 



(2.3b) 



for any pair of momenta q and q' and for any pair of 
functions / and /'. 

The choice f(q,p) = 1 for any pair of momenta q and 
p from M. d defines the momentum representation of the 
local density operator 



p(q) : = j c f (q+p) c(p)- 



Any pair thereof commutes 

[p(q),p(q'))=0. 



(2.4a) 



(2.4b) 



Another choice of the function / is made with the fam- 



ily 



g(q;G):= J e +i ^ q ' p ' G ^ c\q+p)c{p) (2.5a) 



for any pair q and G from M. d where 



§{q,p-G):={ q + G)*p--q*G (2.5b) 



and the ^-product 



a*b = —b*a= a i M^dj 



(2.5c) 



defines a real antisymmetric bilinear form specified by 
the real- valued d x d antisymmetric matrix M'*'. When 
d is even, we assume that is invertible. When d is 

odd, M^*' has at least one vanishing eigenvalue and is 
thus not invertible. Observe that 



p{q) = g(q--q) 

We are going to prove that 



(2.6) 



1. The family g(q;G) labeled by the pair q and G 
from M. d obeys the U(l) algebra (1.1). 



2. In even-dimensional space, the family g(q; G) la- 
beled by the pair q and G from M. d is complete. 



Proof of closure: We define 

r(q,q',p;G,G') := &(q,q' +p;G) 



$(q',pi G')-${q + q',p;G + G') 



in terms of which Eq. (2.3b) can be rewritten as 



[g(q;G),g(q';G')} 



,iT{q,q' , V ;G,G') _ , q ^ g > &nd q ^ Qly 



= +i*(<j+q',p;G+G') -tf_ , „/ 



c\q + q' +p) c(p). 



(2.7a) 



(2.7b) 
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Since 



T(q, q',p; G, G') = (q + \g j * (q' + ±gA - \g * G' = T(q, q'; G, G') 

is independent of p and antisymmetric under q -s-> q' and Gf>G', 

T(q,q';G,G') = -r(q',q;G',G), 



the algebra (2.7b) closes to 

[g(q; G), g(q'; G')] = F(q, q'- G, G') g(q + q',G + G') 

with the structure constant 

F(q, q 1 ; G, G') = e iT( -i<i'-< G - G ') - (q o q' and G -O- G') = 2i sin T(q, q'; G, G'). 
Proof of completeness: Choose any function / : M. d x R d — > C such that the Fourier transform 



(2.8a) 

(2.8b) 
(2.9a) 
(2.9b) 



f(q,p) =-J(<l,P)e iq * P 



3 +i G*p 



f(q,G) 



iq*p 



(2.10) 



is well defined. For the second equality to be true for arbitrary functions f(q, •) 



with the well-defined 



Fourier transform f{q 1 



C, the square matrix M^*> that defines the ^-product must be invertible and thus 



have an even number d of rows (columns). Indeed, the rank of an antisymmetric matrix M^*> is necessarily even. 
Hence, in odd dimensional space, is never invertible as it has at least one vanishing eigenvalue. This means that 

the * Fourier transform J e +lG * p h(G) is at best a function of d— 1 coordinates of p if d is odd. For completeness 

G 

to hold, it is thus necessary that d be even, wh ich w e now assume. A sufficient condition for completeness to h old is 
that the linear space spanned by the operators ( 2.2a ) is limited to the coherent superpositions of the form ( 2.3a) suc h 
that the function /(q, •) : R d — > C has a Fourier transform for any given momentum q. With the help of Eq. ( 2.5b| ), 
we can then write 



f(q,p) = J f{q,Cf)^ mG ' % ^ q ^ 

G 



(2.11) 



In turn and with the help of Eq. (2.3a), we conclude with 



Q f (q) = I f(q,p)£ f (q+p)Hp)= I f(q,G)e [q * G/2 I e^^ c\q + p) c{p) = / f(q, G) e iq * G ^ 2 g G (q)- 



(2-12) 



III. 



THE CASE OF THE LATTICE 



f2 B | +p . Correspondingly, q + p—G q+p € ^ B z- W e shall 
use the notation 



We begin with some notation. Let A be a Bravais 
lattice and A* its dual. Sites in A are denoted by r, sites 
in A* by G. The first Brillouin zone is denoted fi BZ , it 
contains the origin of M. d . We shall decompose M. d into a 



[q + p]bz = q + p G q +p e « 



BZ- 



(3.2) 



Observe that if we hold q € S1 BZ fixed and vary p across 
the Brillouin zone f2 BZ , the unique reciprocal wave vector 



set of shifted Brillouin zones H BZ obtained by translation G q+p G A* such that q + p G 
of BZ by G e A* 



q+p 



€ £! BZ defines 



i d = (J n 



G 

BZ- 



(3.1) 



GeA* 



an implicit function of q that is piecewise constant with 
discontinuous jumps each time q+p crosses the boundary 
separating neighboring Brillouin zones. We define the 
fermionic Fock space J with the help of the algebra 



Sites in f2 BZ are denoted fc, q, and p. If q and p belong 
to the Brillouin zone ri BZ , this might not be the case for 
q + p. There is a unique G q+p € A* such that q + p G 



{ 



} 



(3.3) 
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for any pair k and h! from the Brillouin zone ri BZ and 
any pair G and G' from the dual lattice A*. 

The linear space of fermionic bilinears that we study 
is spanned by the basis 



f •= I st c 



that obeys the algebra 



T ,T , , 



<Ji:Q 2 _ "31.12 I'vll 



(3.4a) 



(3.4b) 



for any quadruple q 1 , q 2 , q[, and q' 2 from the Brillouin 
zone. 

For any q from the Brillouin zone f2 BZ and for any 



function / : BZ x 



"BZ 



define 



9+Wbz P' 



(3.5) 



pe n Bz 

There follows the algebra 

E [/«.ta'+p] BZ 4',p -(<?,/ w,/o 



pen E 
x c [q+9'+p] BZ c p 



(3.6) 



for any pair of momenta q and q 1 from the Brillouin zone 
and for any pair of functions / and /'. 

The choice f q p = 1 for any pair q and p from the 
Brillouin zone defines the momentum representation of 
the local density operator 



E 4* 



"[<J+P] B z 

pe n Bz 



Any pair thereof commutes 



(3.7a) 



(3.7b) 



Another choice of the function / is made with the fam- 



ily 



- E 



[9+PJbz p 



(3.8a) 



p en E 



for any G from the dual lattice and q from the Brillouin 
zone where 



®q.P ■■= ^ [ii + G) *P - (q + P + G) * G q+p+G ] 

(3.8b) 

(*) 

and the d x d matrix ' that defines the *-product 
is antisymmetric, as was the case in the continuum, but 
with the restriction that 



1 

2^ 



to accommodate the d-dimensional Bravais lattice A. 
When d is even, has a nonvanishing determinant 

by assumption. As advertised, the £7(1) algebra 



Qq ) Qq' 



, ( q + G)*(q> + G>^ ~g+G' 



J -q+q' 

(3.10) 

follows for any quadruple q, q', G, and G' . 

To prove completeness, we assume that the dimen- 
sionality d is even for the same reasons as given below 
Eq. ( |2.5c| ). One verifies that 



®q, P = e q , P + G * P 2 /* G + + mod 2tt (3.11a) 



where the function 

2ire qp :=q*p {q + p)*G q+p 
is independent of G, while the function 



2ttQ 



G 



q*G 



is independent of p. We will use the fact that 

®q=0,p = ®q=0 = 



(3.11b) 
(3.11c) 
(3.12) 



in Sec. 1VB and Appendix |A| We define the function 
/ : tt BZ x n BZ — > C by 



f =• f e +ie «." 
We then use the Fourier expansion 



f =• V f 



G +i (G * p— p * G)/(2tz) 

q c 



(3.13) 
(3.14) 



GeA* 



to do the following manipulations 

Jq,p ° Jq,p 



E / 



G +i(G*p-p*G)/(27r) 



G eA* 



- E /■ 

GeA* 

- E ff*- u 



G„+i e qp +i (g * p-p * G)/(27r)+i e°-i e° 



+i * 



(3.15) 



GeA* 



independent of p 



Insertion of Eq. (3.15) into Eq. (3.5) gives 



% = E /^ ie? E 



[q+p] B z p 



GeA* 



■ peiij 



= E f^- 1 



AG 



(3.16) 



GeA* 



where we made use of the definition ( p.8a| ) to reach the 
last equality. Completeness has thus been proved if the 
space of functions / is restricted to those for which the 



G * G' = mod 2tt, VG, G' E A*, (3.9) Fourier transform (pi) exists 
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IV. DISCUSSION 

A pointed out by Murthy and Shankar, the magnetic 
translation algebra is not limited to situations in which 
time-reversal symmetry is broken. From the point of view 
of many-body physics, the generators of the magnetic 
translation algebra can also be thought of as special co- 
herent superpositions of particle-hole excitations above 
the many-body ground state. As such they are always 
present in the many-body Fock space. 

If time-reversal symmetry is either explicitly or spon- 
taneously broken, it is plausible that these excitations 
might be selected by the many-body interactions to play 
an important role at low energies and long distances. 
However, the breaking of time-reversal symmetry alone 
is no guarantee for the FQHE. The selection of a ground 
state supporting the FQHE is a subtle compromise be- 
tween the kinetic energy and the interactions. 

If time-reversal symmetry is neither explicitly nor 
spontaneously broken, it is harder to imagine that these 
excitations are of relevance to the low-energy and long- 
distance properties of interacting electrons. 

We illustrate both observations as follows. We begin 
with the case of interacting electrons in the continuum 
limit without explicit breaking of time-reversal symmetry 
and for which spontaneous symmetry breaking of time- 
reversal symmetry is not anticipated. This situation is 
the one expected if electrons interact through density- 
density interactions. We are going to show how to recover 
the /-sum rule when we choose to represent the many- 
body Hamiltonian in terms of the generators ( [2.5| ) of the 
magnetic translation algebra for any even dimension d of 
space. This exercise serves as a warning for blindly per- 
forming a mean-field approximation of the Ha milt onian 
when represented in terms of the generators (2.5) that 
delivers the FQHE. 

To illustrate the delicate competition between the ki- 
netic energy and the interactions, we consider a band 
insulator to which we add density-density interactions. 
We r epr esent this Hamiltonian in terms of the genera- 
tors ( |3.8| ) for any even dimension d of the Bravais lattice. 
In doing so, we are going to show that normal ordering 
can change the bare band width by a value comparable 
to the characteristic energy for the interactions. Hence, 
if the bare band width is smaller than the characteristic 
energy for the interactions, as is usually believed to be 
necessary to stabilize a FCI, normal ordering is then an 
effect of order one! 



/-sum rule 



U in any dimension d. Pines and Noziercs presented a 
derivation thereof in Ref. ^ that hinges on the fact that 
the operator identity 



Q - 
= — N 



p(q),[H,p(-q)] 
holds for any momentum q € M. d . Here 
N:= J c\p) c(p) 



(4.1b) 



(4.1c) 



is the conserved particle number operator, and the many- 
body Hamiltonian H := H + H v + Hjj is the sum of the 
dispersion 



H Q := J e{p)c\p)c{p), 



(4.1d) 



the one-body potential 



V(q)p(-q), 



(4.1e) 



and the two-body potential 

Hu~ [ V(q) p{q) p(-q). 



(4.1f) 



The only nonvanishing contribution to the nested com- 
mutator in Eq. ( p~Tbl ) arises from the quadratic dispersion 
H in view of the Abelian alge bra ( 2.4b] ). Equation ( 4.1 b| ) 
follows from the algebra (2.1). 

As a sanity check, wc are going to verify Eq. ( 4.1 b| ) 
for any even dimension d with the help of the magnetic 
translation algebra 



[g(q; G), g(q'; G')] = 2i sin T(q, q'; G, G') 
x e(q + q';G + G r 



(4.2) 



where T(q,q';G,G') is defined in Eq. (gj). We shall 
only eva luate the contribution from the quadratic dis- 
persion (4. Id). 

First, we recall that p(q) = g(q; —q) according to 



Eq. (2.6). Second, we expand H in terms of the mag- 
netic translation densities g(q; G), 



e(G) g(q = 0;G) 



(4.3a) 



The /-sum rule holds for electrons with the mass m 
and the quadratic dispersion 

e(p) := ^- (4.1a) 

subjected to any one-body potential V and interacting 
with any translation-invariant density-density interaction 



where 



e(G) 



e(p). 



(4.3b) 



It is with Eq. (4.3b) that we made use of d being even. 
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Third, we make a first use of Eq. (4.2) to evaluate the internal commutator 

'p(q),[H,p(-q)]\ = fi(G)\g(q;-q),[g(0;G),g(-q;q)] 



J s(G) 2i sin T(0, -q; G, q) [g(q; -q), g(-q; G + q)} 



(4.4) 



We make a second use of Eq. (|4. 2f) to evaluate the external commutator 



p(q),[H,p(-q)] 



:(G) 2i sin T(0, -<?; G, q) 2i sin T(q, -q- -q, G + q) g(0; G) 



e(G) 



2i sin 



q* G 



g(0:G). 



(4.5) 



The integral over G can now be performed, 

'p(q),[H,p(-q)]] = fe(G){e +i ^ G + e- i ^ G -2) I e +iG **> c) {p) c(p) 



G 



= / ^ (Ip + if + \P~ if - 2 IP| 2 ) £t (P) KP) 
p 

= to/ at (P) 6 (P)- 



(4.6) 



Equation (4.1b) follows from the definition (4.1c). 



B. Projected Hamiltonians and the importance of 
induced one-body terms 



r 



of sites in A, we can perform the Fourier transformation 
to the band basis in two steps. First, we do the Fourier 
transformation 



We begin with the generic lattice Hamiltonian 



,+ip-r 



H := H n 



(4.7) 



whereby the dimensionality d of the lattice is assumed 
even. Our goal is to understand how does normal order- 
ing of the interaction Hjj changes the band width of the 
kinetic Hamiltonian H . To this end, we need to choose 
the representation in which wc define H and H v . We 
will see that the choice of the representation of H can 
change the effects on H of normal ordering on H v . 
The kinetic Hamiltonian is defined by 



(4.9a) 



in terms of which 



p6Sl E 



p£EQ BZ a, a' 



(4.9b) 



^ E E $r',a> + Hc 

r,r' a, a' 



Second, for any given p from the Brillouin zone, we do 
the unitary transformation 



where the hopping amplitudes 



f ot,ot / i a .a 

r — r 1 1 v' — v 



(4.8a) 



(4.8b) 



lb f =-^C f U a * 



in terms of which 



E^w ( 4 - 10e 



decay exponentially fast with the separation between any 
pair of sites r and r' from the lattice A and we have 
reinstated a finite number of internal degrees of freedom 
labeled by the orbital index a. If N denotes the number 



H - V Vrt f r 

P en BZ a 



(4.10b) 



The algebra (3.3) applies to the band operators labeled 



by the band index a if one multiplies the Kroenecker sym- 
bol S p p , in the Brillouin zone by the Kroenecker symbol 
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8 aa , among the bands. The algebra (3.3) thus endows Normal ordering is the operation by which all creation 



the orbital creation and annihilation operators with the operators are to be moved to the left of the annihilation 
canonical fcrmion algebra. 

The interacting Hamiltonian is defined by 



operators. In the orbital basis, normal ordering results 
in 



with 



r.r' £ A a, a' 

- E EC^^ 

t/Efi BZ a, a' 



rr.a. Yr.a. V^r.a 



(4.11a) 



(4.11b) 



^£7 



(4.13a) 



The one-body Hamiltonian iJ^ , a consequence of the 
fermion algebra, is proportional to the conserved number 
operator, 



the local density at the site r G A and for the orbital a. 
The corresponding Fourier transforms are 



(4.13b) 



r£A a 



JTOt,OL 
1 



N Z^ e ^ ' 



(4.11c) 



For simplicity, we shall focus on orbital independent 
(density-density) interactions, in which case 



where we defined U = U r=0 . The normal-ordered inter- 
action is 



ff V = E Y. U r-r^Ul, a ^r>, a 'i>r, a - (4-13c) 



{/„, Va, a' . 



(4.12) 



The one-body term induced by normal ordering is, in the band basis, 



H 



pGf2 BZ a 



The normal-ordered interaction is, in the band basis, 



(4.14a) 



TjllC 



E 



E E 

<3en BZ pen BzP 'efi BZ 



EEEEC^f 



1+q+p] BZ . a c [-i+p%z^ a ' c p'' b ' C P' 6 ' 



a b a' b' 



v 



a.b:a ,b' 



^E 



u [+<j+p] BZ . a p. & [-q+p'] BZ ' a ' p'. 6 ' ' 



(4.14b) 



In any subspace of the Fock space with a fixed number 
of particles, normal ordering thus produces a rigid shift 
of all single-particle energy eigenvalues of H Q . For any 
band a, the width of the single-particle dispersion e a is 
not affected by the normal ordering — UN. We are going 
to show that this nee ds n ot be true anymore if we first 
project Hamiltonian (4.7) onto the band a and then ex- 
press the resulting projected Hamiltonian in terms of the 



generators (3. 



is 



H a = Hg + Hi 
where the projected kinetic Hamiltonian is 



m= e • 

pesi B 



^p,a ( £ p.a "b ^0 ^p,a 



The projection of Hamiltonian (4.7) onto the band a while the projected interacting Hamiltonian is 



(4.15a) 



(4.15b) 
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E E E 

gGn BZ pen BzP '6f2 BZ 



V a , r , - r 

"q,p,p' °[+q+p] BZ ,a L [-q+p'] BZ ,a V,3 °p,a' 



"q.p.p' — u q "[+q+p] BZ ,a u p,a «[_q+p'] BZ ,a u p\a- 



(4.15c) 



For the p urpo se of representing the projection of 
Hamiltonian (4.7) onto the band a by the magnetic den- 
sity operators (|3.8|), it is necessary to undo the normal 
ordering in Eq. ( |4.15q ). In doing so, a second one-body 
term is produced, 



where the projected kinetic Hamiltonian is 



(4.16a) 



E 



u) 



p - r - 



E 



E 



ya 

q,[-9+p] B z<P 



p,a p,a ! 



(4.16b) 



while the projected interacting Hamiltonian is 



E 



E 

9efi BZ pefi BZ p'en E 



E 



ya »T - -T 

"q,p,p' [+q+p] BZ ,5 u p,a <-[■ 



q+p'] 



,a p ,a- 



(4.16c) 



Observe that had we first represented Eq. (4.11) in 
the band basis, followed by the projection consisting in 
restricting all the band indices to a prior to n ormal- 
ordering, then we would have obtained Eq. (4.16) upon 
normal ordering withou t the second term on the right- 
hand side of Eq. (4.16b). The correct implementation of 
projection is t o norm al order first and then to project, 
leading to Eq. (4.16b). Indeed, the order by which normal 
ordering is followed by restricting all band indices to the 
projected ones corresponds to sandwiching the Hamilto- 



nian by the projection operator onto a subset of bands. 
The reverse order by which the density operators is pro- 
jected onto a subset of bands followed by normal ordering 
corresponds to sandwiching first all density operators by 
the projection operator onto a subset of bands and then 
assembling a Hamiltonian out of these projected density 
operators. As the projection operators do not commute 
with the density operators, the order in which the oper- 
ations of normal ordering and projection arc performed 
matters. 



We can now express the Hamiltonian in terms of the magnetic density operators (the details are provided in 
Appendix |X|) 

H'o a = E \ £ g + US G .o - E ^ ~ h % ) So ' 

G £A* \ q ) (417a) 

H'u a = E % E ff f=? e_i (e? * +e=f ) # *=? . 

gen G,c eA* 
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where 



£ G = 



"q 



J a 



-Ye 
N ^ p 



-i(G*p-p*G)/(2ir) 



N ^ 



+P] B Z' a P' a 



D —i (G * p-p * G) /(2tt) 




(4.17b) 
(4.17c) 
(4.17d) 



with <d q p and defined in Eqs. ( |3.11b ) and ( 3.11c ), respectively. 



Equation ( [4. 17| ) is the main result of Sec. [TV B Ap- 
plied to a Chern insulator to w hich d ensity-density inter- 
actions have been added, Eq. ( 4.17 ) suggests that there 
will always be linear in g G — Q contributions to the Hamil- 
tonian even if the bare band is fiat to begin with, i.e., 
even if e p a = 0. Because of the topological attributes 
of the Bloch spinors as they wrap around the Brillouin 
zone, we expect a nonvanishing h G . (An extreme case of 
a topologically trivial band insulator has Bloch spinors 
that are constant across the Brillouin zone, in which case 
only h G=0 7^ and the additional one-body contribu- 
tion is just proportional to the total particle number. 
This would also be the case in the context of the quan- 
tum Hall effect.) This effect on the bare dispersion is 
controlled by the bare interaction U q . Hence, it could be 
as large as the effects of the density-density interaction. 

It is far from evident that a FCI is selected by the 
competition between the one-body and two body terms 
in Eq. (4.17) since they are both controlled by one char- 



acteristic energy scale in the limit of a flat bare band 
width. On the other hand, if a ground state supporting 
an FCI is selected for some range of parameters, then the 
effective quantum field theory describing the low-energy 
and long-distance properties of this phase should belong 
to one of the universality class associated to the FQHE. 
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Appendix A: Details for reaching Eq. (4.17) 



Equipped with Eq. (4.16), we are in position to take 
advantage of the fact that, for any even dimension d, 
the magnetic density operators (3.8) form a complete ba- 



sis for the charge neutral fermion bilinears made out of 
the band creation and annihilation operators. Define the 



functions /* : f2 BZ x f2 BZ — > C, with i = 1, 2, 3, 4 by 
f 1 ■= /) r 

J q.p ■ q,0 p,a ! 
fq,p := 3q,0 U- 
f, 



lq',p' 



q,p 



~ fiq',0 ^p',p> 

[q+p] B z> s p. 5 ' 



(Ala) 



(Notice that p is a parameter in the d efin ition of / 3 .) In 
terms of the operators defined in Eq. ( |3.5D , the projected 
kinetic Hamiltonian is 



n/ 1 



E Kt- 



qen E 



p,q'GO BZ 



while the projected interacting Hamiltonian is 



(Alb) 



(Ale) 



The bare kinetic energ y is the first term on the right- 
hand side of Eq. ( Alb ). The correction to the bare ki- 
netic energy from standard normal o rderi ng is the second 
term on the right-hand side of Eq. ( Alb ). The function 
f 2 is a delta function with respect to its first argument 
and constant with respect to its second argument. This 
correction does not change the band width of the bare 
dispersion. The last correction to the bare kinetic energy 
is the third term on the right-hand side of Eq. (Alb). It 



is controlled by the interaction U q dressed by the Bloch 
functions that diagonalizc the bare kinetic energy, 



ya 

q,[-<?+p] B z>P 



-<3+P] B z> a U P' a 



(A2) 



It can change the bare band width by an amount of the 
order of the ch aract eristic energy of the interaction. 

At last, Eq. ( |3 .13 ) dictates that we define the functions 
f : n BZ x n BZ -> C, with i = 1, 2, 3, 4 by 



Jqp 



f 
■I a 



(A3) 



where the argument of the multiplicative exponential fac- 
tor on the right-hand side was defined in Eq. (3.11c). We 
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observe that any of the three fun ction s f 1 , f 2 , and / 3 is 
proportional to 8 q so that Eq. ( 3.12 ) teaches us that 



fq.p : ~ ^q,0 /o,pJ * — 1; 2, 3. 



(A4) 



In even dimensional spac e, we can safely use the *-Fourier 
expansion [see Eq. (3.14)] 



flp = S q ,o E fo G e+^ G *P-P* G y^\ i = 1,2,3, 
GeA* 

14 . \p ?4;G +i(G*p-p*G)/(27r) 
■'Q.P Jq <=■ ) 



E 



G eA* 



E 



-G +i(G*p-p*G)/(27i-) 



GgA* 



(A5) 



to compute the Fourier coefficients f^ ,G with i = 1,2,3, 



fq' G , and h G . Application of Eq. ( 3.1 6| ) then delivers 
the d esire d representation of the projection of Hamilto- 
nian (|4.7|) o nto the band a by the magnetic density op- 
erators (|3.8|) , 



3£A* \ q J 



G £A* \ q 

//;- E E .r ; /V ; < i(0? « +e O ?■ 

qen BZ G,c ga* 

(A6) 



Equation (4.17) follows. 
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